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Abstract: 

We consider the polarization of A + A baryons produced in polarized Deep Inelastic 
Scattering at leading order, with various spin configurations: longitudinally po- 
larized leptons and unpolarized nucleon; unpolarized leptons and longitudinally or 
transversely polarized nucleons; longitudinally polarized leptons and nucleons. We 
show how the different results in the different cases are related to different aspects of 
the elementary dynamics and to the spin properties of the distribution and fragmen- 
tation functions and show how a combined analysis might give useful information. 
We give numerical results according to several sets of polarized fragmentation func- 
tions recently proposed. 



1. Introduction 



A baryons produced in liigti energy interactions and resulting from quark frag- 
mentation allow a unique test of spin transfer from partons to hadrons: the A po- 
larization is easily measurable by looking at the angular distribution of the A ^ pvr 
decay (in the A helicity rest frame) and the fragmenting parton polarization is de- 
termined by the elementary Standard Model interactions, provided one knows the 
initial parton spin state. In this respect A's produced in lepton induced processes 
are particularly interesting and indeed several papers on this subject have recently 
been published or submitted to e-Print archives 

We perform here a detailed analysis of the polarization of A's produced in po- 
larized DIS; a general discussion of the helicity density matrix of hadrons produced 
in polarized lepton-nucleon interactions, at leading order, can be found in Refs. |1lO| 
and [|I|. It can be written as: 

rlrr^,S+N,S^t+h+X 1 

(h) = E ^ X (1) 

'^^h^^h'- ' dxdydz ,V V levrxs ^ ' 

pX,\ P\^K fll^^^^ ^A„A,;A„A, ^A* A;;A^,A; ^ ^"^J^z) , 

where x and y are the usual DIS variables, x = Q'^/2p-q, y = Q'^/xs and, neglecting 
hadron masses, z = Ph'P/p-Q, where p, q and ph are, respectively, the nucleon, virtual 
photon and final hadron four-momenta, p^'* is the helicity density matrix of the 
initial lepton with spin s, fq/N{x) is the number density of unpolarized quarks q with 
momentum fraction x inside an unpolarized nucleon and the helicity density 

matrix of quark q inside the polarized nucleon with spin S. The ^^ .^^^ ^ 's 
are the helicity amplitudes for the elementary process iq — > £q. The final lepton 
spin is not observed and helicity conservation of perturbative QCD and QED has 
already been taken into account in the above equation: as a consequence only the 
diagonal elements of p^'** contribute to p{h) and non diagonal elements, present in 

A A' 

case of transversely polarized leptons, do not contribute. D^' y {z) is a generalized 

h ' h 



fragmentation function related to the usual unpolarized fragmentation function 
Dh/q{z), i.e. the density number of hadrons h resulting from the fragmentation of 
an unpolarized quark q and carrying a fraction z of its momentum, by 

DH/q{z) = ^ E ^a':;a:(^) = I E , (2) 

A,,A^ A,,A^ 

where -D^' a'' (^) = hx (^) is a polarized fragmentation function, i.e. the density 

h h Q 

number of hadrons h with helicity A^^ resulting from the fragmentation of a quark 
q with helicity A^. Angular momentum conservation and collinear configurations 
imply for the generalized fragmentation functions ||10|| : 



d',1''^=0 when A,-A;^A,-A;. (3) 
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The elementary amplitudes M"^ are normalized so that the elementary unpolarized 
cross-section, £q — in, is: 



2/ 



(4) 



Notice that there are only two independent helicity amplitudes: 
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(6) 
(7) 

dy dy 167ra;s ^'^ Q2 y (S) 

Finally, the cross-section appearing in the l.h.s. of Eq. (|l|), which gives the 
correct normalization to p{h), Trp = 1, can be written, using Eq. 



so that 



with 



and 
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{y)\'D,/g{z) . (9) 



Eqs. hold within QCD factorization theorem at leading twist and leading 

order in the coupling constants; the intrinsic k± of the partons have been integrated 
over and collinear configurations dominate both the distribution and the fragmen- 
tation functions. For simplicity of notations we have not indicated the scale 
dependences in / and D. 

We shall use such equations for spin 1/2 A baryons produced starting from 
several particular initial spin configurations; we will discuss how the measurable 
components of the A polarization vector depend on different combinations of distri- 
bution functions, elementary dynamics and fragmentation functions: each of these 
terms predominantly depends on a single variable, respectively x, y and z, and a 
careful analysis of different situations can yield precious information. Although A 
production in polarized DIS has been recently discussed in several papers, most of 
them only consider some initial spin configurations and specific models for fragmen- 
tation functions. Our analysis is more comprehensive and somewhat more general, 
emphasizing the physical meaning of possible measurements, and allowing also to 
obtain some general relationships between different polarization values. 
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2. Polarization vector of spin 1/2 baryons 

We fix our spin notations in tlie I — p center of mass frame: tlie lepton moves 
witli four-momentum I along tlie z-axis and the proton moves with four-momentum 
p in the opposite direction; we choose xz as the lepton-hadron production plane, 
with the y-axis parallel to / x p^. We denote by 5*^ the (longitudinal) nucleon spin 
oriented along the z-axis, by Ss the (sideway) spin oriented along the x-axis and by 
Sm the (normal) spin oriented along the y-axis. Notice that -1-5'^ corresponds to a 
— helicity proton and — S'j;, to a + helicity one; we will only consider longitudinally 
polarized leptons with spins is/, which correspond respectively to ± helicities. 

From Eqs. ®)-® one obtains the explicit expression for the components of 
the helicity density matrix of a spin 1/2 baryon. These are related to the three 
components of the baryon polarization vector, as measured in its helicity rest frame 
by Pi{B) = Tr[((T*p(-B)] {i = x, y, z); details can be found in Ref. [jlO|- We denote by 
{s,S) [or {h,H)] the (lepton, nucleon) spins [or helicities], stands for unpolarized 
particle; one finds, for a spin 1/2 hadron B: 



Y.qe\ Ag(x) ADBiq{z) 
Ege2g(x) DBiq{z) 



(10) 



E,(l/2)g[rfa++-rfa+-] AD, 
E,e2g(x) ADBiq{z) 



Eg q da'i UB/q 



J2q el q{x) DB/q{z 
P^''''^\B-x,y,z) = -P^'''''^\B-x,y,z) 



AlM (11) 
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P^^' ^^>{B;x,y,z) = P^+'+>{B; x,y, z) = —— ! / 

Eq[q+ d(T++ + g_ da+ \ DB/q 

Eqcl [q{x)ALL{y) + Aq{x)] ADB/q{z) 



Ege2 [q{x) + ALL{y)Aq{x)] DB/q\ 



z 



(13) 



Pi--^)(E;.,.,.) ^ Pi^-) (B; X, y, z) ^ ^j}]' ^^^'^ 

T.q[q~ da++ + g+ da+ \ DB/q 

Y.qel [q{x)ALL{y) - Aq{x)] ADB/q{z) 
EqCl [q{x) - iiL(2/)Ag(x)] DB/q{z) 

where d&q stands for ddq/dy. 
Some comments are in order. 



(14) 
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The above results are known []T], |n[]; we have rederived and grouped them 
here for convenience and further discussion. Notice that all other spin config- 
urations - at leading order - either give the same results or no polarization. 
We remind that qx = fq^/p+ is the number density of quarks with helicity A 
inside a + helicity proton, q{x) = + q-{x) and Aq{x) = q+{x) — q-{x). 

Similarly AD^/q = Dsj^/q^ — DB_/q+', A^g and AxDs/q are respectively the 
analogue of Aq and ADs/q for transverse spins. 

The longitudinal B polarization induced by a longitudinal nucleon polariza- 
tion, Eq. (0), does not depend on the elementary dynamics, but only on the 
quark spin distribution and fragmentation properties. Neglecting the QCD 
Q^-evolution, P^^'^^l) (Joes not depend on the DIS variable y, but only on x 
and z. 

The longitudinal B polarization resulting from the scattering of longitudinally 
polarized leptons off unpolarized nucleons depends on the unpolarized distri- 
bution functions, the polarized fragmentation functions and the elementary 
dynamics, through the double spin asymmetry for the £q iq process [see 
Eqs. (i-®]: 

Notice that An grows with y from (at y = 0) to 1 (at y = 1), so that P^^l,o) 
is an increasing function of y, starting from at ?/ = 0. 

The transverse B polarization induced by a transverse nucleon polarization, 
Eq. (p!2D, depends on the quark transverse spin distribution and fragmentation 
properties and on the elementary dynamics, through the double transverse spin 
asymmetry for the iq iq process: 

DnnW -- 



da'i'-^i' - d&^'^'^'i' 2(1 - y) 



da'i 1 + (1 - yf 



(16) 



where t= Sn and [= —Sn- 

Contrary to An, Dnn decreases with y, with Dnn = 1 at ?/ = and Dnn = 
at y = 1; thus, pj'^'-^^) is a decreasing function of y, reaching at ?/ = 1. An 
experimental confirmation of the opposite y-dependences of P^^lA) g^y^d pio,SN) 
would supply a new, subtle and important test of the factorization scheme of 
Eq. (|I]). Dnn is the transverse polarization of the final quark generated by 
an initial transversely polarized (f) quark in the iq iq process, and it is 
usually referred to as the depolarization factor. 
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• When both the lepton and the nucleon are longitudinally polarized the B re- 
sulting polarization depends on yet a different combination of polarized quark 
distribution functions, fragmentation functions and elementary dynamics. It is 
then clear why a combined study of P{B) in different cases could yield unique 
information. 

• Finally, we notice that, in general, P^'^'^^ = — p/ ^\ 



3. Polarization vector of A baryons 

We now consider the particular case of A baryons and discuss possible ways 
of extracting information from combined measurements of -Pi(A); as we said the A 
polarization vector can be measured by looking at the proton angular distribution as 
resulting from A np decay in the A helicity rest frame, that is the frame obtained 
by rotating the i — p cm. frame around the |/-axis so that the new z/j-axis is parallel 
to the A direction, and then boosting along with the same speed as the A: 

W{9p,(j)p) = — [1 + a;(P2 cos + Pj, sin 6'p cos 0p + Pj^ sin 6'p sin 0p)] 

= ^l^ + ^P-P] (17) 

where a = 0.642 ± 0.013. 

We follow Ref. 0] and assume for the unpolarized fragmentation functions: 

Da/u = Dk/d = Da/s = Da/u = Da/S. = Da/s = Da/q (18) 

where A means A° + A". 

The heavy quark and gluon unpolarized fragmentation functions play a negligible 
role for z > 0.3 and we neglect them here: we have actually checked that our 
results, when comparable, are almost indistinguishable from those of Ref. P] where 
also heavy quark and gluon contributions to the unpolarized cross-sections are taken 
into account. 

Similarly, we follow Ref. for the polarized fragmentation functions: 

ADA/uiz, Ql) = AD^/diz, Ql) = ADa/,{z, QI) . (19) 

Eqs. (|19D holds also for light antiquarks and it remains valid through QCD Q^- 
evolution; heavy quark contributions are neglected. 

Using Eq. (||) and (|l|) into Eqs. ©-([ill) and (|ll)-(|l) gives: 



Q{x) DA/g{z) 

P. [A,x,y,z} - ^^^^ ^^^^^^^ AllIv) (21) 
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where 



Pi--)(A;x,,,.) = Q'i-)^^^(y)^^Q'(-) ^^A/;(f (22) 

Qix) + AQ{x) ALLiy) D^,^{z) 

' ^ ' ^ Q{x)-AQix)ALLiy) D^,^{z) ^ ^ 



Q = 4{u + u) + (d + d) + (s + s) (24) 

AQ = 4(Am + Am) + (Ad + Ad) + (As + As) (25) 

Q' = [A{u + u) + {d + d)]N^ + {s + s) (26) 

AQ' = [A{Au + Au) + {Ad + Ad)] + {As + As) . (27) 



If one assumes the same relation (|T9D to hold also for transversely polarized quark 
fragmentation functions, then Eq. (|T^) yields: 

p{o,5^)rA - ^ „ A - ^tQ'{x) AtDa/s{z) ^ , , 

with 

AtQ' = [4(Atm + Atu) + {Ard + Ard)] + {Ats + Ats) . (29) 

Eqs. (^0])-(|^) hold under assumptions (|l^) and (|19D alone, independently of the 
actual value of A^^; as the polarized and unpolarized distribution functions are well 
known and several sets are available in the literature, we can exploit Eqs. (^OD and 
(plD to obtain: 

AT — - 

U All - (Af//f/) P] 

and 

ADa/s{z) _ Q 1 {AU/U) Pi+-°) - All 
Da/,{z) S All i^U/U) - {AS/S) 
where we have defined 



^ _ S All - {AS/S) Pj+^o) 



(31) 



U = A{u+u) + {d+d) S = s + s AU = A{Au + Au) + {Ad+Ad) AS = As + As. 

(32) 



Eqs. ( p2D and ( 123|) can then be used to predict - within the general assumptions 
(|18D and (|T9|) - the following interesting relations between polarization observables: 

p(+,0) ^ p(0,+) 

p(+,±) ^ -^^ =^ (33^) 
1±(AQ/Q)iii 

and 

1-[{AQ/Q)All? ^ ' 
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Notice that in the small y region, due to the elementary dynamics, see Eq. ([T3|), 
one has: 

Pi+'°) ~ Pi+'±) ~ ±Pi°'+) < 1) . (35) 

In the large y region instead, again from Eq. (|15D and from the fact that large 
y implies small x, where Q{x) ^ AQ(a;) and Q'{x) ^ A(5'(x), we expect: 

p{+,+) ^ p{+,-) ^ pj+,0) ~ 1) . (36) 

We conclude this Section by reminding that we have derived our results for 
A = A° + A°, which allows assumption (|TB|) concerning q and q fragmentation func- 
tions. However, Eqs. (p0| ) -(p3| ) and (^) hold identical also for single A° production, 
provided one neglects the fragmentation function of a g into A°, i.e. one neglects all 
q terms in Eqs. (p4D-(p7D and (pQ]). Anyway, the production of A°, in processes 
is strongly suppressed by the hmited amount of initial g, unless one considers very 
small X values. 

4. Numerical estimates 

We give now some numerical estimates of Eqs. (^D|)-(^^ and (^); we use the 
sets of unpolarized and polarized fragmentation functions introduced and discussed 
by the authors of Ref. 0]: together with Eqs. (|18|), we use the expression for the 
unpolarized fragmentation functions they obtained by fitting e^e~ KX data. At 
initial Ql = 0.23 (GeV/c)^ scale one has, from a leading order (LO) analysis |^: 

Dj,/,{z,Ql) = 0.63z'-''il-zY-''. (37) 

The polarized fragmentation functions are assumed to be of the initial form 
dlD, with: 

ADa/,{z, Ql) = z''Da/,{z, Ql) . (38) 

Leading order QCD evolution is consistently taken into account in our numerical 
computations.^ Next to leading order contributions to the A polarization have been 
shown to be tiny and we neglect them. 

The parameter Nu defined in Eq. ([T9|), has been chosen according to three 
different scenarios typical of a wide range of plausible models, and the corresponding 
remaining parameter a of Eq. ( pS]) has been fixed by fitting the few LEP data on 
A polarization, with the results 0: 

1) Nu = 0, a = 0.62. This scenario corresponds to SU{6) non relativistic 
quark model, according to which the whole A spin is carried by the s quark. 

2) Nu = —0.2, q; = 0.27. Such value of Nu is suggested in Ref. ||12[|, based on 
a SU (3) flavour symmetry analysis and on data on the first moment of gf. 

^We are very grateful to D. de Florian, M. Stratmann and W. Vogelsang for providing us with 
their FORTRAN package for unpolarized and polarized fragmentation functions. 
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3) Nu = 1, q: = 1.66. A scenario in which, contrary to the non relativistic 
models, all light quarks contribute equally to the A polarization. 

The unpolarized and polarized distribution functions are taken respectively from 
Refs. [jl3[ and ||14[. We have explicitely checked that a different choice of unpolarized 
and polarized distribution functions, like those of Refs. |jl5| and [Q, can change 



significantly the numerical values of -Pj(A) (but not their qualitative behaviour) only 
for X > 0.3. 

A computation of Pj'^''^'^-', Eq. ([28|), requires the knowledge of the quark trans- 
versely polarized distributions, AtQ or hiq, and of the transversely polarized frag- 
mentation functions, AtD, which are not known. In order to give an estimate we 
fix Axq for u and d quarks by saturating the Soffer's bound [|17| (assuming the same 
signs for Arq and Aq): 

Atu = ^{u + Au) Ard = {d + Ad) . (39) 

All other transverse distributions (A^^ and Aj-s) are neglected here and we also 
assume At-Da/s = AD^/g. 

A sample of typical results is presented in Figs. 1-4, for HERMES kinematics, 
s = 52.4 (GeV)2 and > 1 {GeV/cf. 

• In Fig. la we plot Pj°'+) - at fixed x = 0.1 and 2; = 0.5 values - as a function 
of y, for each of the three scenarios; P^^'~^\ Eq. (pOD , can depend on y only via 
the Q^-evolution and indeed the three curves show an almost flat behaviour. 
The three scenarios yield quite different results. The minimum value of y is 
given by ymin = Q'LiJ^^s) - 0-19. 

The same plot for pj+'*^) is presented in Fig. lb; the y-dependence is essentially 
due to the factor Ai^i^ in Eq. (pi]); scenario 3, which assumes Q'{x) = Q{x), 
together with Eq. (|38D in which we neglect the mild Q^-evolution (taken into 
account in our numerical computations), gives a particularly simple result: 

p(+,0)^^1.66^^^_ (40) 

• Pi"*"'^-* and P^^~\ again as functions of y at fixed a; = 0.1 and z = 0.5 values, 
are shown in Figs. 2a and 2b respectively. At large y ^ 1 values Eq. (^) is 
satisfied. The small y — >■ behaviour cannot be seen with x = 0.1; in Fig. 2c 
and 2d we plot respectively P^^'~^^ and P^~^~\ changing the x value to x = 0.3 
and keeping z = 0.5. The allowed minimum value of y is now 0.06 and we 
have checked that Eq. (pS]) is indeed obeyed (by comparing with Pi°'+) as a 
function of y at x = 0.3 and z = 0.5, not shown in Fig. 1); the change in sign 
of Pj^'"-* is particularly interesting. 

• In Figs. 3a-3d we show respectively Pj^'^-*, Pz^'^\ ^z^'^^ -^i^' ^ fixed 
values of = 1.7 (GeV/c)^ and z = 0.5, as functions of x. At fixed Q^, y 
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decreases with increasing x - and viceversa - and this explains why relations 
(^) and (^) hold respectively when x ^ 1 and x — >■ 0. Once more, the three 
different scenarios give very different results. Using different sets of polarized 
distribution functions, like |]15| or [|T6[ instead of |T^, gives almost identical 



results for x < 0.3 and larger (in magnitude) results for x > 0.3 (of course, 
Pj"*"'") is not affected at all by a change in Ag). 

In Figs. 4a we plot P^^'^'^\ Eq. (pSD , at fixed x = 0.1 and z = 0.5 values, as 
a function of y for all scenarios; the y dependence is almost entirely given by 
Dnn, Eq. (0), and indeed pJo,5iv) _^ g -^j^en y ^ 1. Notice the opposite y 
behaviour of Pj+'°) and P^(0''S'jv) ^j^g ^i^e opposite behavior of All and D^n- 

In Fig. 4b P^°'^^^ is plotted function of x at fixed = 1.7 (GeV/c)^ and 
z = 0.5. 

In general one finds very small values of Pi (A) in scenario 1, negative values 
in scenario 2 and positive ones in scenario 3. This can easily be understood 
by the different values of in the three scenarios, which assign respectively 
zero, negative, and positive contributions to u and d quarks, which dominate 
in the proton. Experimental measurements can easily discriminate between 
them. 



Recently HERMES Collaboration [18| have published a single experimental 



measurement of P^'^'^y All, as a function of z and this seems to favour the 
scenario 1 prediction of Ref. 0, although errors and uncertainties are still 
too large to draw any reliable conclusions. Similarly, two values of P^~^'^\z) 



published by the E665 Collaboration [O] still have much too large errors. 



5. Conclusions 

The study of the angular distribution of the A decay allows a simple and 

direct measurement of the components of the A polarization vector. For A's pro- 
duced in the current fragmentation region in DIS processes, the component of the 
polarization vector are related to spin properties of the quark inside the nucleon, to 
spin properties of the quark hadronization, and to spin dynamics of the elementary 
interactions. All this information, concerning quark distribution functions, quark 
fragmentation functions and spin properties of elementary dynamics are essentially 
factorized and separated as depending on three different variables, respectively x, z 
and y. The Q^-evolution and dependence of distribution and fragmentation func- 
tions somewhat mix the three variables, but smoothly, keeping separated the main 
properties of each of the different aspects of the process. Moreover, such depen- 
dence is perturbatively well known and under control. 

We have discussed all different polarization states of baryons, obtainable in the 
fragmentation of a quark in DIS with polarized initial leptons and nucleons, Eqs. 
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(|T^-(|n|), showing how they can reveal different quark features, weighted and shaped 
by elementary dynamics. 

Adopting a simplifying - although rather general and representative of many 
possible choices - assumption about the quark fragmentation functions into a A Q], 
we are able to extract from measurements further information on the quark frag- 
mentation process, Eqs. ( pOD and (0), and to predict relations among polarization 
states induced by different initial spin configurations, Eqs. (|33D and (0). 

Numerical estimates are given in Figs. 1-4, according to three largely different 
scenarios for fragmentation functions; each scenario has physical motivations 
and yields qualitatively different results: compatible with zero, large and negative, 
large and positive. Such results are stable against different choices of the polarized 
and unpolarized distribution functions, so that experimental information should 
immediately allow to draw clear conclusions and to learn about quark fragmentation 
properties. 

The elementary dynamics fixes the small or large x or y behaviour of some of 
the polarization components; although expected, such behaviours should indeed be 
checked, as an independent and non trivial test of the QCD factorization scheme 
of Eq. (HI); such a scheme has been widely used and tested for the computation of 
semi-inclusive unpolarized cross-sections, but not for more subtle spin observables. 

We think that our comparative and comprehensive discussion of all possible A 
polarization measurements in polarized DIS is useful and can lead to a new and 
clear strategy which allows to obtain novel information. 
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Figure 1: Pj'^'+^ and Pj"'"'^) as a function of y at fixed values olx — 0.1 and z — 0.5, 
for the three different scenarios. 
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Figure 2: Pi+'+) and Pi+'-) function of y at fixed values oi z = 0.5 for the 
three different scenarios. The upper plots correspond to fixed a; = 0.1 and the lower 
plots to X — 0.3. 
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Figure 3: Pi°'+), Pi+'°), P^+'+^ and P^+-^ as a function of x at fixed values of 
= 1.7 (GeV/c)^ and z = 0.5, for the three different scenarios. 
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Figure 4: (a), P^^'^^) function of y at fixed x = 0.1 and 2; = 0.5; 
(b), P^(0'^^) as a function of x at fixed = 1.7 (GeV/c)^ and z = 0.5. 
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